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We consider the effect of the curvature in fermionic dispersion on the observable properties of 
Luttinger liquid (LL). We use the bosonization technique where the curvature is irrelevant per- 
turbation, describing the decay of LL bosons (plasmon modes). When possible, we establish the 
correspondence between the bosonization and the fermionic approach. We analyze modifications 
in density correlation functions due to curvature at finite temperatures, T. The most important 
application of our approach is the analysis of the Coulomb drag by small momentum transfer be- 
tween two LL, which is only possible due to curvature. Analyzing the a.c. transconductivity in 
the one- dimensional drag setup, we confirm the results by Pustilnik et al. for T— dependence of 
drag resistivity, R12 ~ at high and R12 ~ at low temperatures. The bosonization allows for 
treating both intra- and inter-wire electron-electron interactions in all orders, and we calculate exact 
prefactors in low-T drag regime. The crossover temperature between the two regimes is Ti ~ -EfA, 
with A relative difference in plasmon velocities. We show that A 7^ even for identical wires, due 
to lifting of degeneracy by interwire interaction, U12, leading to crossover from R12 ~ U?2T'^ to 

R12 ~rVt/i2 at r~(7i2. 

PACS numbers: 73.63.Nm, 71.10.Pm, 72.15.Nj, 73.23.Ad 



I. INTRODUCTION 

The effects of the curvature of fermionic dispersion for 
the observables in strongly-correlated one-dimensional 
(ID) fermionic systems were discussed by several groups 
recently jii^i^i'^i^i^i One physical effect, which probes both 
the curvature and interactions is the Coulomb drag effect. 
In a typical experimental setup, this effect is observed as 
non-zero drag resistivity, R12 = V2/Ii, with a d.c. cur- 
rent Ii flowing in the " active" wire and the voltage bias 
V2 per unit length apphed to the second "passive" wire 
in order to assure J2 = 0, see review articles Ref. Q. 

It is generally agreed that the principal source of the 
drag effect is the particle-hole asymmetry of the elec- 
tronic system. This was established in various works con- 
sidering electrons in higher spatial dimensions, D > 1, 
both in the absence of magnetic field^ and in the pres- 
ence of it (seei^ and references therein). It was shown 
that the leading contribution to the drag effect is ob- 
tained in the second order of the interwire interaction 
U12 and is schematically depicted by the Feynman graph 
in Fig. [2l This contribution corresponds to the virtual 
processes away from the Fermi surface, and is ultimately 
determined by the curvature of the electronic dispersion. 

One naturally would expect that the curvature should 
be also responsible for the drag effect in one spatial di- 
mension (ID). However, the first study of the drag effect 
in ID was devoted to a different mechanism^ii specific 
for ID and requiring some additional conditions, we re- 
turn to it below. The study of the drag effect in ID due 
to curvature of electronic dispersion was initiated only 
recently, by Pustilnik et ali The reason is that the drag 
effect depends both on the curvature and the fermionic 
interaction, whereas it is known that the interaction is 
very important in ID and ultimately leads to the notion 



of Luttinger liquid. The possibility to find a complete 
solution, i.e. Luttinger liquid, for the interacting system, 
lies in the crucial simplification of the theory - the lin- 
earization of the fermionic dispersion around the Fermi 
energy. In most cases it is an innocent procedure, and 
does not infiuence the final result. The Coulomb drag 
effect is one example, where this usual theoretical trick 
with linearization leads to the immediate disappearance 
of the observable quantity in question. In their work^ 
Pustilnik et al. (see also^) hence used the fermionic for- 
malism, and found the drag effect as the function of tem- 
perature, for different regimes. Restricting themselves to 
the lowest order of interwire, as well as intrawire inter- 
actions, they showed that R12 (x for equal wires and 
i?i2 oc for non-equal wires. In this work we check their 
results by another method and extend their treatment 
beyond the lowest order of perturbation theory (PT) in 
fermionic interactions. 

To this end, we use the traditional tool for studying the 
effects of interaction in ID, known as bosonization tech- 
nique, which is devised for effective resummation of the 
appearing series in PT.-^^ Employing the bosonization we 
take into account the inter- and intrawire interaction in 
all orders, whereas the curvature of fermionic dispersion 
is treated as perturbation. We show that it suffices to 
consider fourth order in curvature when calculating low- 
temperature drag effect. At higher temperatures we per- 
form the resummation of most singular contributions in 
PT in curvature, incorporating free-fermion results into 
bosonization language. In this sense, our method is com- 
plementary to the analysis ini where the curvature was 
treated exactly and the interactions — in the lowest nec- 
essary order. 

The issue of the curvature of fermionic dispersion in 
bosonization technique was rarely discussed previously. 
Therefore, although we are ultimately interested in the 
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drag effect, we develop a systematic approach to Lut- 
tinger liquid with curvature. It was shown long ago ^^i^^ 
that the curvature corresponds to the cubic terms in 
bosonic densities (plasmon modes) of so-called right- and 
left-moving fermions near two Fermi points ±kp, respec- 
tively. The appearance of these cubic terms corresponds 
to the decay of bosons (plasmons), and was previously 
derived in two ways via operator identities. We suggest 
yet another explanation for the appearance of the cubic 
terms in the theory, an explanation referring to Green's 
function formalism and more pertinent to our discussion 
of Coulomb drag. Besides that we discuss the modifi- 
cation of the current operator, dynamic density correla- 
tions in isolated wires, and the light-cone singularities at 
finite temperatures T. Keeping T finite is important for 
our calculation, as the drag effect is only due to inelastic 
processes of fermionic interaction, and hence vanishes at 
T = 0. We use the Kubo formalism, which technically 
delivers the optical transconductivity, aij{uj), in one di- 
mensional system of two wires, and eventually extract 
the d.c. drag resistivity, a~j^{0). Given our assumption 
of clean Luttinger liquids, it is a delicate procedure and 
we discuss it in detail. 

The plan for the paper is as follows. We set up the 
problem and discuss the identification of fermionic cur- 
vature with decay of plasmon modes in SecHT] The form 
of density correlation functions is discussed in Sec. Illli 
we show the correspondence between fermionic approach 
and bosonization here. The bosonization expression for 
current and conductivity for one wire is discussed in Sec. 
IIVI The transconductivity matrix for two wires is con- 
sidered in Sec. |Vl and the drag coefficient is obtained for 
relatively simple case of non-equal wires. The analysis of 
nearly-equal wires is performed separately, in most tech- 
nically involved Sec. IVI[ the final results for R12 in low 
temperature regime are discussed here. We present our 
summary and conclusions in Sec. IVIll 



II. SETTING UP THE PROBLEM 
A. Hamiltonian 

We consider spinless (spin-polarized) electrons in two 
wires, with the forward-scattering short-range interac- 
tion both inside and between the wires. The electrons 
in ith wire (z — 1,2) are conventionally subdivided into 
right- and left-going species, centered around correspond- 
ing Fermi momenta, kpi ^ Trn^, with the fermionic den- 
sity Ui. We go beyond the the linearized spectrum ap- 
proximation, keeping quadratic terms in dispersion. 

As usual, we decompose fermionic operator into 
"slowly oscillating" chiral components 

^(x) = e'^^-^Vflla;) + e-^'^-^M^), (1) 
so that the smooth part of the fluctuating fermionic den- 



sity is 

p{x) = ^lj(x)^ij(x) + ^pI{x)Mx). (2) 
Our Hamiltonian is given by four terms : 

n = J dx{Hi+H2+ H12 + Hcur), 

with the linearized interacting fermions in individual 
wires 

+U,p%x)l2, (3) 
forward-scattering interaction between wires 

Hl2 = Ui2Pl{x)p2{x), (4) 

and the terms describing the curvature of fermionic dis- 
persion 

Hear = X! {i'Rji'idxf'^Rj+i'Lji^dxfi'Lj) ■ 

. (5) 

The bosonization approach to the above system is the 
application of the important representation of ID chiral 
fermions as 

V'ij, = -i=e^^«^ Vi,— ^^e-'^-^ (6) 
\/2^ \/2^ 

with A the ultraviolet cutoff. Here he bosonic fields 
4'R(L),j = 4'3 T ^ji with primary field (pj and its canon- 
ically conjugated momentum H^ = ir^^dxdj satisfying 
['/'j(2;),n,(y)] = iSjiS{x-y). 

In this bosonization notation^^^, the Hamiltonian of the 
system is rewritten as follows : 

n = J dxiHi+H2+Hi2+Heur), (7) 

Hi2 = — ^9^0ii9x02, (9) 

Hcur - (3(9,0i)2 + (a,0i)2) (1 ^ 2). 

DTrmi 

Here and below I use the shorthand notation dx(f>i = 
d(l)i{x)/dx, etc. The electronic density operator in ([7]) is 
given by pi — dx4>i/T^- We follow the notation by Hal- 
dane, with Vij = vpi = kpi/mi, ViN = vpi + Ui/n. In 
the absence of i?i2, Hcur, the dispersion of (plasmon) 
excitations is Si{q) — Vi\q\, where vf — VijViPf. The 
strength of the intrawire interactions may be encoded in 
the Luttinger parameters, Ki — y^Vij/viN ; it is conve- 
nient also to define the dimensionless interwire interac- 
tion, U = Tr~'^Ul2/y/VlNV2N- 

It is seen that the Hamiltonian ([7|) contains an exactly 
solvable quadratic-in-bosons part Hi + H2 + H12 and the 
curvature acts as the interaction. 



B. Curvature as interaction 

In bosoiiization, the normally ordered parts of chiral 
fermionic densities become 
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(10) 



with pj — Rj + Lj 



^dx4>j- The main statement 



of (Abelian) bosonization, dating back to Tomonagar^ 
is that the free chiral fermions with linear dispersion 
are equivalent to quadratic form in chiral densities : 
—itp]idxipfi = ttR^. Tomonaga established this relation, 
studying the equation of motion for the chiral density R 
with the fermionic Hamiltonian with linear dispersion. 
This observation was later corroborated by Schick,— 
who found that similar consideration with the quadratic 
fermionic dispersion yields another important relation : 
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(11) 



Another derivation of (|lip was shown in the paper by 
HaldanCfi^ who started with the fermion representation 
dni). Note that using Eq. one arrives also to the term 
in the Hamiltonian, which constitutes a full derivative, 
~ d'^(j)fidx<pR ~ dx{R^), which does not contribute to the 
dynamics. Confining oneself with the quadratic terms in 
the fermionic dispersion, one obtains nothing beyond the 
terms R^ + and R^ + in the bosonization language. 
The stability issue for the cubic action is briefly discussed 
in the next section. 

Keeping in mind a problem of drag effect, let us provide 
yet another argument verifying the form of the curvature 
term in bosonization. It is well known that the multi-tail 
fermionic loop diagrams in one spatial dimension are ex- 
actly zero for the linear dispersion law. This statement 
is known as the higher-loop cancellation theorem after 
the work by Dzyaloshinskii and Larkinji^ One can show 
however, that the presence of the curvature in fermionic 
dispersion leads to generally non-zero value for the pro- 
cesses given by triangle diagram. The expression for this 
diagram, characterized by three external vertices, can be 
easily calculated at T = as 
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A(w,q) 



{uJiq2 ~ L02qiY 



1 



In 



u) — vpq - 



- (91 9293/2^)2 ' 
q-^/2m 



2tt lu — vpq — q'^ /2m 



(12) 
(13) 



with (73 = qi + q2, LU3 = uji + 102- One can check that 
T vanishes in the limit m — s- cx). Expanding (|12p to the 
leading order of m~^, one obtains the structure 



1 



92 



93 



277771 Vpqi — OJi fF92 — OJ2 fF93 — OJ3 ' 



(14) 



which corresponds to three density correlators for lin- 
ear spectrum, {2Tr)~^q/ (vpq — (^), attached to above ver- 
tex, 27r^/(3m); factor 3! = 6 comes from symmetrization. 
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FIG. 1: Feynman diagrams representing non-linear response 
of fermionic system. Two three-tail fermion loops correspond 
to a boson decay process, see the text for discussion. 





FIG. 2: First Coulomb drag diagrams in the fermionic for- 
malism and their bosonization counterpart, see the text for 
discussion. 



Schematically, it is shown as Feyman diagrams in Fig. [U 
where fermion Green's functions are shown by lines with 
arrows, and double lines represent bosonic Green's func- 
tions ((2T|) . Note that the leading-order drag diagram is 
given by the "two triangles" fermionic diagrams in Fig. 
[21 with wavy lines standing for Ui2^^ Comparing it to 
Fig. [1] we see that the drag coefficient should correspond 
to "two stars" bosonic diagram in the r.h.s. of Fig. [2l 
Further observing that interwire interaction, U12, is ex- 
actly included into the Luttinger liquid formalism, Eq. 
([7]) , we eventually arrive to the diagrams depicted in Fig. 
m below. Note that passing from Fig. [2] to Fig. [4] corre- 
sponds to the screening of interwire interaction, see e.g. 
Kamenev and Oreg^^ The RPA resummation in case of 
LL gives an exact result; however, the simplest bosonic 
diagram in Fig. [2] corresponds to the optical transcon- 
ductivity, rather than to d.c. drag effect, as will be shown 
below. 

Before considering a general case of two different wires, 
let us discuss the modification of density correlations in 
one wire with the curvature term. 



III. DENSITY CORRELATIONS 



Let us first discuss the density correlation function for 
one wire, with and without the Luttinger-liquid type in- 
teraction. Such analysis was previously done in'^''*, we 
extend it for the case T 7^ 0. The bosonized Hamiltonian 
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reads as 



W = y dxiHo + Hcur), 



Hq = — (vj{dx0f + VN{dx(, 
Ztt 

U 



(15) 
(16) 



7rwF(i?2 + + -(i? + L)2, (17) 



1 



dx(p {z{djf + {dx4>f) , (18) 

(19) 



67rm 

- 2^ 
3m 

The commutation rules for the chiral densities are 



+ m 



[R{x),R{y)] = -[L{x),L{y)]^—dJ{x-y).{20) 

ZTT 

Our choice of the interaction in the form U{R + L)"^ cor- 
responds to the case 174 = g2 in more general form of the 
Hamiltonian 34 (i?^ + i^) + 2g2RL. 



A. Free fermions 

In the absence of the mixing interaction term URL in 
(fT7|) . each chiral component is considered separately. In 
the absence of the curvature, 1/m = 0, the bare boson 
retarded right (left) Green's functions is independent of 
temperature and given by 



(0) 

9r{l) 



1 



27r qvF =F (w + iO) ' 



(21) 



with (5— function spectral weight. 

Importantly, in this free case we know the exact ex- 
pression for correlation function of the right (left) mov- 
ing densities. The retarded (right) density propagator 
x{i-^,q) can be calculated at non-zero temperature as the 
fermionic loop diagram. 



X 



dk tanh(efe/2r) - tanh(efc+q/2r) 



47r 



(22) 



with £k = vpk + /2m. Neglecting the terms of order 
of T/{mvp), we get : 



Im X 



sinh(u;/2T) 



4g cosh ( '^^^^g-^^ ) cosh (gj 



(n+i; 



4T 



(23) 



AlTWa 



^ ' 2 AttT 



AttT 

{ll> ~ qvF)/wq, Wq^q^/{2m), 



(24) 
(25) 



with ip{x) digamma function. Restoring the real part of 
X we used the rule 



; tanh ; 



iIj ( i — 

TT^ \ 2 TT 



for the function analytic in upper semiplane of z. In the 
limit of zero temperature we obtain 



Im X 



X 



?e(i 

2q ^ 



AnWa 



In- 



qvp 



Wo 



(26) 
(27) 



J ~ qvp 

with the step-function 8(a;) = 1 at x > 0. 

Let us now discuss how these exact expressions re- 
late to the bosonization theory. Formally, the curvature 
term H^ur has a scaling dimension 3, which is strongly 
irrelevant. The common wisdom then would prescribe 
to treat it in the lowest order of perturbation theory 
(PT) in the infrared limit, g, w 0. Comparing ([27|) 
and the bosonization expression ((2T|) . one can see that 
the (5— function is a good approximation for the spectral 
weight in this limit. However, the lowest PT contribu- 
tion, with two bosons in intermediate state, leads to the 
self-energy part of the form 

" "2 ^ (2^tIvfY 



E = 27r 



12to2 



(28) 



qvp — uj — iQ 

Note, that in a quite unusual wayi^ the spectral weight 
Im S(a;, q) is located exactly at the "light cone" of bare 
spectrum, uj — qvp — 0. The Dyson equation, = 



9r 



S, reads as 



= 2n 



qvp — uj 

q 



{2ttT/vp) 



21 



12m2 



qvF 



(29) 



It means that instead of previous single pole ui = qvF^ we 
have two poles separated by a distance ~ m~"'^g maxfg, T] . 

Now we compare (p9)) with (|24|) . by expanding x~^ in 
powers of . We find that ([^ is reproduced by ([M]) in 
the order m^'^ . It suggests that bosonization produces a 
correct asymptotic series in m~^, but at the energies close 
to "light-cone" condition, an ever increasing number of 
PT corrections should be taken into account, to arrive at 
the exact answer 

Interestingly, the bosonization predicts also that the 
fermionic interaction of the form gi{R^ + L'^), i.e. with 
g2 — 0, results only in the shift of the Fermi velocity, and 
therefore the spectral weight x has the previous form ([M]) 
except for the change, vf ^ vf + 9^121:. This statement 
should be contrasted with a recent work^ where a change 
in X in this case was estimated by resumming the RPA 
series with the 54 term. However, in the presence of 
curvature, the RPA series does not exhaust all possible 
processes, as the multi-tail fermionic loops do not vanish, 
Eq. Hence the result of simple RPA resummation is 

different from the answer by bosonization, which implies 
the multi-tail fermionic processes, particularly, the "drag 
diagrams" , Fig. [21 



B. Curvature in Luttinger liquid 

Let us discuss now the modification of the theory 
in the presence of the mixing interaction ^ g^RL in 
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(fT7| . In this case right and left chiral densities mix, but 
the quadratic Hamiltonian which is the Luttinger model 
is still solvableji^ The corresponding u — v Bogoliubov 
transformation diagonalizes the Hamiltonian in terms of 
new chiral densities, R, L, which obey the same com- 
mutation relations pop . Performing this transformation, 
one obtains 



Ho 


= TTV 


{k^ + 


Hcur 


_ Ctl 

3 


[R^ + 






3 + K^ 








m 


2VK ' 



(30) 



-fiR'L + Rp), (31) 



a2 



7r2 (^2 



m 



K 



(32) 



with — vf{vf + U/i:), = (1 + U/nvp)'^- One can 
see from (|3T|) that apart from some rescaling of the old 
mass term, W + L^, the effect of interaction gives rise 
to a new type of the curvature vertex, which mixes new 
right and left densities. 

The density correlation function {pp)uq is determined 
through the new variables i?, L as 



{pp)^^ = K{{R + L){R + L)) 



Note that RL mixing, which was removed at the level of 
quadratic action, reappears due to the curvature. The 
mixed terms (RL) are non zero in case of both U ^ 
and ^ 0, and it necessitates the consideration of 

the (retarded) matrix Green function 



pressions: 

Arr 
All 
Im Arl 



1 g(g2 + (27rr/z;)2) 
487r'^ qv ~ u! ~ iO 

1 q{q^ + {2nT/v)^) 
487r3 

2 



qv + UJ + iO 



2 



? 9 r 



(47ru)3 

oj — vq 



coth 



UJ + vq 
AT 



coth 



4T 



(37) 
(38) 

(39) 



Note that the anomalous (appearing for K ^ 1 only) 
contribution Arl is ultra-violet divergent, which will be 
reflected in the appearance of UV cutoff A ~ Ep in real 
part of Arl- Seeking a function Arl, analytical in the 
upper semiplane of w, one can restore Re Arl using the 
rule 



i coth f3u! 



.f3uj 



1 

i(3uj 



Such a recipe of analytical continuation delivers Re Arl 
up to a constant ^ ln(T/A) multiplied by (w^ — v^q^), as 
will be seen shortly. 

In the limit of T = 0, the expression for Arl is sim- 
plified 



Im Arl = TT- 

Re Arl = - 



v^q'^ 



(47r-(;)3 



Qiuj"^ ~v^q^)signuj, (40) 

2 



(47rv)' 



q^v^ 



A2 



.(41) 



G{uj q) = ( ^^^9 



It is clear that the PT for bosons in the case K ^ \ 
contains processes absent in the free fermion case K = 1. 

Let us calculate the boson self-energy in the second 
order of PT. Denoting three types of boson loops as 
Arr, Arl, All (Fig-)> one obtains 



^RR — 

^LL - 
^RL — 



2alARR + -alALL + (^Arl, (33) 



2a\ALL + ^alARR + alARL, (34) 
T.LR^ aia2{ARR + All) + ct^ARL- (35) 



Off the light cones, oj ^ ±vq, one has Im Arr = 
Im All — and all components of the imaginary part of 
the self energy are given by the quantity 



Importantly, this part of S vanishes at the light cone for 
T = and remains finite (and purely imaginary) for T ^ 
0, in contrast to the singular contributions Arr, All- 
Notice the appearance of UV cutoff A in Eq. (|4ip , which 
replaces temperature in this limit. 

Neglecting for a moment the singular (purely real) 
parts Arr, All, we write the Dyson equation in the form 



G- 



9r^ ~ ^off 



(42) 



with g^i^ given by (|2T|) with the change vp ^ v. Some 
calculation shows then that the Green's function defined 
for the initial fields (j>, 9 is given by 



\{de,d(j)), {80,89} 



7rg2 / (^/q 

D^t \uj/q, K-\v~^off/^)) ' 

(43) 



with 



So// = a\ARL. 



(36) 



2„,2 



Det = q^v 



q^v 



We list here the form of the appearing (retarded) ex- 



and 8 standing for spatial derivative in (I43p . 
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As a result, the density correlations are given by 

Kvq^ /t: 



{PP)u 



n 

VF 



w2 - q^vYjoff/n' 



In the limit T = we have 

Z-^ = l-7(9/fcF)'ln 



2 2 

V q 



A2 



7 



64 



(44) 
(45) 

(46) 
(47) 



and kp = mvp — Kmv. The obtained logarithmic cor- 
rection to the residue Z should be regarded as a first 
term in a series in m~^. This correction is increasingly 
important at the energies close to light-cone condition, 
uj v\q\, and corresponds to the first term in the ex- 
pansion of power-law singularity at cu ~^ v\q\, which is 
discussed in^. 

Away from the light cone, when the real part of loga- 
rithm in Z can be neglected, we obtain the contribution 
to the spectral weight at w > 0: 



Im {pp)ujq 



7 '^-F'? 



Y^Q{lo' - v'q'), (48) 
qZyZ 



which shows a tail at higher energies, > \qv\. This 
result was obtained earlier iit^'^. 

Finite temperatures, T > \lu\, \qv\, lead to a following 
modification : 



2 

q V , 



(49) 



and from ([44]) it appears that the bosonic excitations are 
characterized by damping 



±vq + 2i'KjT{q/kp)^ 



(50) 



However, this "damping" should not be viewed as the 
characteristic linewidth, because it is smaller than the 
previous estimate for the width of the spectral width of 
density correlator in the noniteracting case, \uj — vq\ ^ 
Tq/kp- This wider spectral width comes from the omit- 
ted terms Ahji^All which intervene close to the light 
cone condition lo — ±vq. The finite-T "damping" due to 
interaction manifests itself in the enhanced amplitude of 
the tails in the density correlator, and instead of ([^5]) we 
obtain 



Im (pp) 



7 



vpq^ 



Luq 



2kj, {ui' 
+ coth 



LO — vq 



coth 



vq 



AT 



7 



AT 
Avpq^LoT 



kp (w^ — (p-v^Y 



(51) 
(52) 



which is parametrically larger than (|48)) at \jJ — vq\ ^ 
LO T ^ and is non-zero also at |a;| < \qv\. It should be 
stressed again, that Eqs. ((48|) . ([52|l are applicable only 
for \lo — wg| > fc^^l^l max(T, qv). 



IV. A.C. CONDUCTIVITY OF SINGLE WIRE 

Let us discuss the coupling to electromagnetic field in 
bosonization and its modification due to curvature. We 
start from general expression for the kinetic term in the 
Hamiltonian, 

■Hk^n = 2^^t (-iV - eAf^, - V, (53) 

with ti — c — \ and p chemical potential. Using the 
above representation V-" (a;) ^ gjfcF2;+»0H_|_g-ifcFa;-»0i,^ ^nd 
denoting 

^R^'kpX ^ ^L = kpX + (l)L 

one can represent Ti.kin in the presence of eA as 



-N - 1- 

ZTT 



{d.,<^R - eAf 



1 

'2^ 



6 m 

{d,<^L + eAf 



6m 



pdx^R 



(54) 



Slightly digressing here, we note that in this notation 
kpx is clearly a zero (classical) mode of the bosonic fields 
^R,L which is obtained from ([51)1 by variation over these 
fields. Note that the zero mode kpx corresponds to a 
local minimum of the action corresponding to (j54p . while 
generally speaking the action is unstable with respect 
to choice dx^R,L — ^ — oo. The latter unphysical choice 
of " classical" vacuum can be ruled out by imposing the 
hard-core condition, dx^R,L > for all x. It means that 
the local fluctuations of the density should not lead to 
negative values of total electronic density, dx<t>R > —kp 
etc. 

Note, that we can obtain ([M]) also by including vec- 
tor potential in the U{\) phase according to ipix) ~ 

^t-S>n-^eJdxA _^ ^-^'^I^~^e j' dx A ^ rpj^g ^^^^^^ ^^^^ -^^ 

generated by two terms in (|54p cancel each other. In more 
general case of non-parabolic band dispersion, e{k), the 
Pcicrls substitution, k ^ k — eA, leads to higher powers 
of eA in the gradient expansion of e{k — eA). The recipe 
for tackling this case was outlined long ago by Luttinger 
and Kohn^® and amounts to retaining only eA and {eAf 
terms in such expansion, while discarding higher terms 
as unphysical. 

Upon thermodynamic averaging of TLkim the diamag- 
netic term, ~ {eAf, gives the usual value, kp/{'Km) = 
Kv/tt. This is the zero- mode contribution, whereas 
terms linear in dx4> yield zero. The term linear in vector 
potential has a structure 



dx 



^Lf - {kp + dx<jyRf] 



TT '^eA{vpdxO + m ^dx0dx(i 



FIG. 3: Two boson loop diagrams formally contributing to 
a.c. conductivity. The labels 6(j>, etc., at boson lines refer to 
the boson fields as disussed in text. 



FIG. 4: Two boson diagrams for a.c. conductivity, obtained 
with the use of continuity relation and equivalent to Fig. |3] 



The last expression means that the electric current is 
given by the expression 

3 = -eTT-\vFdJ + m~^djd^(l)). (55) 

This expression is confirmed by the immediate appli- 
cation of the continuity equation, dxj = —dtp — 
ieTT~^[dx4>jT^]i a-iid using the bosonized Hamiltonian 
P^ . The first term in (|55p is usual for LL, the second 
term appears due to fluctuation of the Fermi velocity 
when changing the density. 

The Kubo formula expresses the conductivity via the 
retarded Green's function of the uniform current, 

( N {ij)q=a,u - {jj)q=0,oj=0 

crfcjj = . (56) 

lUJ 

The subtraction of the diamagnetic term {jj)q=o,uj=o 
should be performed when using direct representation 
([55| . If one calculates first the density correlator and 
uses the continuity equation afterwards, 

then the diamagnetic contribution is subtracted auto- 
matically. 

Using the above result, Eq. we have 

a(uj) = — , Re ct(w) = vfS(lu), (57) 
—muj 

i.e. the conductivity of one wire is not affected by the 
temperature, interaction and curvature, in accordance 
with Giamarchi and Millis.iS, It should be stressed that 
this conclusion assumes our above choice 52 = 54 = 
i.e. physical charge-charge interaction ; we do not con- 
sider models with 172 ^ 54 which would correspond to 
current-current interaction and to modifications of (j57p . 

We also comment here on the equivalence of two ways 
of calculation of <j{uj), via Eq. (j55p and via continu- 
ity equation. Using (|55p for (jj) correlation function, 
one obtains a bosonic propagator {dxO, dxO), and bosonic 
loops of the form {dx9dx4', dxOdx(t>) ■ From Eq. (HS)) 
the average {dxO, dxO) is vanishing in the limit q — > 0, 
and subtracting the diamagnetic contribution wf/tt, one 
arrives at ([57)1 . The bosonic loop terms give zero, as will 
be seen shortly. 

The mixed terms in {jj) correlator, m~^{dx0,dx0dx4') 
are zero in all orders of perturbation theory. It is trivial 
to see in lowest order of m~^, because the mixed aver- 
age contains three bosons. In the next order, one has a 



diagram with a bosonic loop and the Greens's function 
attached to it. This Green's function bears the zero ex- 
ternal momentum and finite uj and hence vanishes. The 
same argument about zero contribution to conductivity 
can be applied to any diagram containing external lines 
carrying zero momentum, particularly in the two-wire 
situation. Hence, apart from the (S— function contribu- 
tion produced by conventional bosonization, we have to 
discuss only loop diagrams in Fig. [3] and their dressing 
by the interaction, m~^. 



V. A.C. CONDUCTIVITY, TWO WIRES 

For the system described by the Hamiltonian ([7]), we 
wish to compute the transresistivity, i.e. the voltage drop 
in one wire, caused by the current in another wire. Tech- 
nically, we use Kubo formula and calculate the conduc- 
tivity, not resistivity. We consider the drag effect which 
is caused purely by electronic interaction, and in the ab- 
sence of disorder in wires. In this situation, one expects 
a S— function contribution to the conductivity, which 
should be separated from the proper drag contribution, 
seen in the limit lu 0. Let us first discuss this small-cj 
limit. 



A. d.c. limit 

A general form of transconductivity in the limit w — !■ 
can be found, when starting with the transresistivity ma- 
trix R{uj) = [d-{uj)]^^. Exactly at = the voltage drop 
in both wires should be absent when the drift charge ve- 
locities, TTji/kpi, in both wires coincide. In order to see it 
suffices to make a Galilean transformation to the coordi- 
nate frame moving with the drift velocity.^ We combine 
this observation with the third Newton law, equalling the 
forces acting on both wires, i.e. Eikpi = —E2kp2, and 
represent the complex- valued impedance as a series in uj 

i?(^) = .ro (58) 

-iujA + uP'Ai + . . . , 

with the transresistivity given by i?i2 = t^tqI (}ip\kp2)- 

In the absence of interwire interaction U\2 in ([7]), 
Eq. ([nil) shows that ro = , A2 = and A = 
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TT diag[vpl,Vp2]- We make a natural conjecture, that 
rQ,A,A2 allow Taylor expansions in u. It follows then, 
that for Lu < rovjkp^ the first two terms in ([58|) are of the 
order of rgA;^^, and the third term is ^ (rofc_F)~^0(u), 
i.e. contains an additional small factor O{uro). Therefore 
one can invert R{uj) at small u, retaining only first two 
terms in ([58]) and taking the above non-perturbed value 
of i. 

The resulting complex-valued conductivity is obtained 

as 



a{uj) = 0-1(0;) -f (72 (w) 

,.-1^-1 

CTl (w) 



-1-0 - icj \^ r7ij^^m2""^, Vp\vF2'mi'^ 



0-2 {uj) 



1^-1 



-2 



rod — V~"^i ^"^2 



(59) 



(60) 



seek a lowest order contribution to Re <t{uj) behaving 
as w"^. Such contribution should arise in the higher- 
than-second order of m^^, and the next important term, 
uj^^, is found in even higher order of PT. Our approach 
in determination of as the coefficient before the uj~'^ 
term in conductivity is thus similar to memory function 
formalism. 

Below we prove that apart from the singular cu = 
contribution, the matrix ^{uj) is given by 



(64) 



This degenerate matrix cannot be inverted and one needs 
the complete representation ([S^]) in order to obtain Ri2. 
From now on we discuss the real part of the drag con- 
ductivity ad{(jj)- 



Interestingly, (|59p provides the universal, i.e., indepen- 
dent of details of fermionic interaction, relation between 
the drag resistivity and the drag conductivity in small-w 
limit. 



i?i2Re [CT12] U^o = {c + c 



(61) 



with = mikpi/ {•m2kp2). The r.h.s. of (j6ip is always 
less than 1/4, which value is achieved for equal wires. 

In the continuum media, the integral of cr(w) over all 
frequencies should not depend on the interactions in the 
system, which statement is known as the optical sum rule 
(see, e.g.j^ and references therein). From (|59p we restore 
this statement : 



Re / dujd{ijj) 

-O(rod) 



VFl, 
0, Vf2 



(62) 



which means that the spectral weight of the part (72 (w) 
is large ^ vp and is concentrated at the lowest ener- 
gies. Regrettably, this is the worst accessible region for 
bosonization analysis, as will be seen shortly. 
We can cast (1591) in the form 



+ 



VFl, 
0, VF2 

u){u> + irod) \—m 



(63) 



and identify the first term here with the bare bosons con- 
tribution ([57]) . then the second term is entirely the effect 
of the curvature. At larger co, we have the expansion 



Re a 



Tom 



rod 



The drag effect, r^d ^ 0, is absent without curvature, 
and bosonization treats the curvature as perturbation. 
It means, particularly, that calculating the d.c. drag re- 
sistivity starting from Kubo formula, (j56p . one should 



B. Optical Transconductivity 

For the Hamiltonian ([7]) , the bare Green's function de- 
fined for the vector <I> = [dxdi, dxipi, dxd2,dx<p2) is found 
as 



G = TT 



'vij, u/q, 0, 

Lo/q, WiAT, 0, Uu/tt 

0, 0, 'i;2j, ^/g 

0, Ui2/7T, U!/q, V2N 



(65) 



It shows that the excitation spectrum has two branches 
E-j-q = v+\q\ and e-q = V-\q\ with the sound (plasmon) 
velocities 



4 = l{vl+vl)±^\{v!-vir + u^vlvl (66) 

we assume w+ > w_ > below. We also focus on the 
case of almost identical wires, vi V2 — v and small u, 
when ~ «_ ~ u and 



A = {v+ - v^)/v^ ~ y/{vi - V2y/v'^ u2 < 1. (67) 

Let us discuss the first diagrams appearing in bosoniza- 
tion analysis of transconductivities. As shown above, the 
linear-in-bosons components of the currents (j55|) . d9i, 
produce the ^— function contribution to the conductiv- 
ity. This is the first term in Eq. ([63| . 

The proper transconductivity part, (7d(w) in (j64p is 
generated by the second term of the current operator in 
([55)) . Schematically, it is given by two diagrams, depicted 
in Fig. [3| These bosonization diagrams correspond to 
the usual fermionic ones. Fig. [2l after we identify Fig. [2| 
with Fig.[3|in the approach using the continuity relation, 
and notice that the external tail Green's function become 
constant in the limit q 0, u ^ 0. 

After a lengthy calculation, sketched in Appendix, we 
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obtain 



„wsinh(w/2T) vlvlv~\z\v+ + v^y^ 

11 — 

+ ^ (68) 



/ UV1V2 \^ v+ + S-y+w^ 



(69) 



First observation about the obtained optical conductivity 
is that (Td(w) non-zero solely due to the interwire inter- 
action u. Secondly, at lu > T the second term in (I68p is 
exponentially small and the first one produces linear-in-o; 
optical conductivity, 



(70) 



Unusually, there is no upper boundary for this regime in 
the theory, and Read{oj) oc |a;| until |u;| < Ep. At the 
same time, the integral contribution of higher frequencies 
into the optical sum is small, / duj u(ui) ~ v^vp, i-e. only 



a fraction of the above value 

In the case of identical wires, vi = V2 = v, coupled by 
a weak interaction, one has u <C 1 and the new velocities 
v± — v^l ± u. The optical conductivity takes then the 
form 



o'd(w) 



u'^uJC0ih[u}/4T) 
128^3 



wsinh(w/2T) 
8u3u3 sinh^(w/2wr)' 



(71) 



One sees the appearance of a smaller energy scale uT <C 
T in this case and parametric enhancement Ud{Q) ^ T/u^ 
we return to this point below. 

Overall, the function ([55)1 is a smooth function of w, 
not exhibiting the internal crossover at w ~ rgd ^ m^^ 



as suggested by ([63|. The region of applicability of the 
expression can be estimated by using the "broad- 
ened" propagators ([M]) instead of 5— function- like (pijl . 
during calculation of crd(ti'). This estimate shows that 
the finite bosons' linewidth intervenes the formula ((68)) 
at w < Wo = T^/mu^, in accordance with^^. Using 
this estimate as a crossover energy from to the 

Lorentzian form (|63p . we can write crd(wo) — ''o/'-^o and 
thus find up to a numerical coefficient. Remarkably, 
this simple method gives a correct value of ro at lowest T, 
but it wrongly estimates a crossover to high-temperature 
regime, discussed below. We see that the estimate of tq 
requires going beyond the lowest order of curvature and 
we suggest a way to calculate it in the next subsection. 



C. Alternative representation of (Jdiuj) 

In this subsection we derive another representation for 
CTd(w), which explicitly recovers the matrix structure Ij64p 
and shows the overall prefactor v? in Udioj)- For infinites- 
imal u this representation yields the well-known formula 
for the drag coefficient. 



We rewrite the Kubo formula as 



cr(w) 



(72) 



and use the equation of motion dtj = i\H,j] and Eq. 
(j55p . Some lengthy, but straightforward calculation ver- 
ifies that the commutation with the terms Hi,H2, Hcur 
produces full derivatives with respect to a;, which vanish 
for the uniform currents. The only surviving terms come 
from H12 in ([7]) and are 

dtji{x) = Ui2m^^pidxP2, (73) 
dtj2{x) = Ui2m2^P2dxPi ^ -Ui2'm~^pidxP2, (74) 

modulo full derivatives. It is immediately seen that the 
structure (j64p is reproduced, and the scalar drag conduc- 
tivity is given by 



[/2 

O'd(w) = —^{pidxP2, PldxP2)q=0,t 

Comparing to we obtain at a; > r^d 



r/2 



■Im {pldxP2: PldxP2) 



(75) 



(76) 



For the infinitesimal U12 the cross-terms, ~ (pi, dxP2), 
in the average (|75p can be neglected, and a comparison 
with Eq. at w ^ gives 



q^dqduj Im xi{q,uj)lv[Y X2{q,^) 



(2^)2 



2Tsinh^(w/2T) 



(77) 



which corresponds to previous findings. ^'^ It is also worth 
noting that our Eq. (j76p for transresistivity exactly corre- 
sponds to Eq. (18) in^S,. We should stress, however, that 
the exact formula for the case of clean Luttinger liquids 
is given by Eq. ([75]) . and (|76|) should be regarded as the 
high-energy asymptote of the expression 



D. Non-equal wires 

We consider now the case of two wires, which are cou- 
pled by infinitesimal U12 and are characterized by non- 
equal plasmon velocites. In this case we can use Eq. 
([77)1 and the results of the Section IIII Bl In the absence 
of curvature and intrawire interaction, the formula (I44p 
with Z = \ gives tq = in (|77p . The first non- vanishing 
contribution is given by (j52p . and we obtain 



U?2VpiVF2 



q^dq 



lenTivl - v^) J sinh i^^^ sinh ^^^^^ 



72 



Vikp2 sinh(uiq/2r) 



+ (1^2) 



(78) 



In the important case of almost identical wires, when 
A = (vi — V2)/v ^ 1, we can omit the wire index in the 
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prefactors, expand sinh((7Ji 
get 



V2)q/AT) ~ Avq/iT and 



(79) 



with 7 given by (|47p and we used dx / sinh^ x — 
7r^/15. Comparing ([79]) to the result by Pustilnik et air- 
we see that we obtained the same dependence but 
parametrically larger prefactor, {vi — W2)~^, rather than 

Note that one would reproduce the prefactor obtained 
ini, if one neglected the T— dependence of density prop- 
agator and substituted in ([77]) the expression (|48p for 
Im xi(g,a;) calculated at T = 0, instead of ([5^ . We 
see that finite T parametrically enhances the "tail" of 
density correlator, and it has consequences for the drag 
effect. The importance of calculation of intermediate po- 
larization operators Im XiilT'-^) in (|77p at finite T was 
also emphasized 

Checking the applicability of ([5^ . one can conclude 
that the derivation of ((79|) assumed a condition vi — V2^ 
T/kp. In other words, the expression ([7^ should be 
regarded as low-temperature regime of drag resistivity, 
taking place at 

T<Ti = kpii'i - V2) = EpA. 

Above this temperature scale, the " core" spectral weights 
of Xiil^) a-iid x'iil'^) overlap and one can estimate the 
drag coefficient, by taking the expression for the "core" 
of the density propagator for free fermions, Eq. to 
obtain 



T^F 



2T 



(80) 



with 



F{z) 



2x dx{x coth X — z coth z\ 
Q sinh(a; -|- z) sinh(a:; — z) 
= 1, z^O, 
~ 4z3e"273, z > 1. 



The estimate shows an activational behavior of the 
drag for free fermions, similarly to Ref. [l|. 

The obtained estimates for the drag effect at lower 
and higher temperatures, however, raise a few questions. 
First, the Eqs. ([50]) and ([7^ apparently describe differ- 
ent contributions to the drag, they do not match at the 
suggested crossover temperature T ^ Ti, and Eq. l|75|) is 
explicitly non-zero only due to intrawire interactions and 
this feature is absent in ([50)1 . 

Second, in case of initially equal wires, the plasnion 
velocities are split due to interwire interaction, Eq. (|67p . 
It means that U12 cannot be regarded as infinitesimal in 
([75)1 and hence Eq. ([75]) is, generally speaking, invalid. 
In this case a naive substitution A ^ U12/V from ((67|) 
results in ro, Eq. ([75)1 independent of interwire coupling. 

To clarify these issues, we separately consider the case 
of nearly equal wires below. 



VI. DRAG FOR NEARLY EQUAL WIRES IN 
BOSONIZATION 



A. Identical wires, Hamiltonian 

When the wires are identical, it is convenient to intro- 
duce the symmetrized combinations 



± , 



±( 



(81) 



A technical convenience of \/2 prefactor stems from the 
same commutation relations for the symmetrized oper- 
ators, however the zero modes are different, dx<j)^^^ = 

V2kFj dx'j)'"^^ = 0. We will call ^+ and (j)- pseudocharge 
and pseudospin fields, by obvious analogy with bosoniza- 
tion for one wire in the spinful caseJ^ 

In terms of densities and momentum operators, 

p± = 7r-ia,0±, n± = n-^d,e±, (82) 

the Hamiltonian (O is rewritten as follows 

n = J dx{H++H^+ H,.,r), (83) 
H± - ^ivjnl + VN{l±u)pl), (84) 



Hrii.r ' ^ 



{p\ + 3n^ + an^ + ipl) + 6n+p_n_ 



6m\/2 



It is seen that the change in the zero mode of 4>+ is com- 
pensated by the enhancement of its "mass" . The plasmon 
velocities and Luttinger parameters are given by 

= {l±u)vNV,j, Kg^ = [1±u)vn/v,j. (85) 

According to our definitions 

K^^ = 1 + {tivf)-\Ui±Ui2), 

so that = 1 for equal intrawire and interwire inter- 
actions. 

We diagonalize the Hamiltonian ([83]) by transforma- 
tion 



p± = k]/^p±, n± = K^'''U±. (86) 

In terms of new left and right movers, 2R± = p± + H±, 
2L± = p± — H-t, there are no mixed terms R±L± in 7i, 
and propagators {R±L±) vanish. 
The curvature term becomes 

Hcur = (-^+ + -^+) + 'Y (^+'^+ + -^+-^+ 



1/2, 



+az{R+ + L+)R^L^, 



(87) 
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FIG. 5: Skeleton diagrams depicting the bosonic processes 
contributing to drag effect at lowest temperatures. 




FIG. 6: RPA processes, which do not contribute to drag 
conductivity. 



with 



ai 



"3(4) 



as 



Kl 



0L2 




m 



V2 



1 2 

\- K ± 

mV2 VX- K 



iV2 



1 



in the limit of vanishing interwire interaction we have 
— K- and 2ai = and 02 — — as] this case 
returns us back to Eq. (|3ip. In the absence of interaction, 
K± = 1 and a2 = 0:4 — — 0. Note that the vertex 
ai, describing the decay within the same mode, is always 
non-zero. 

The drag current operator in (I75p is 



pidxP2 = P-dxP+ = [K+K-Y^'^p-dxP 



(89) 



modulo full derivatives. We seek the contribution to ro 



in the next order in 



which can be obtained by con- 



sidering two skeleton graphs, depicted in Fig. [5l The 
diagrams of the RPA type. Fig. [6l vanish, because in our 
case the intermediate propagator, connecting two boson 
loops, carries the zero external momentum and is zero at 
any external w 7^ 0, see (|2T|). 

Given ten possible decay processes (|57)) . one has to 
analyze ~ 10^ diagrams, corresponding to Fig. [51 It is 
further complicated by the necessity to perform the cal- 
culation at finite T and make an analytical continuation 
to real frequences w — > 0. Fortunately, there is a method 
which allows the essential simplification of this task. 



B. Calculation by "Unitarity condition" method 

At finite T, the analytical continuation of the diagrams 
and the essential reduction of their number under con- 
sideration are simultaneously achieved in the "Unitarity 
condition" method.— Since this method for T > is 
largely unknown to the wider audience, we sketch the 
main points of it below. Readers not interested in the 
details of calculation may skip to Eq. (|98p below. 



FIG. 7: Four cross sections corresponding to two above skele- 
ton graphs. 



We define a cross section of the diagram as a line, bi- 
secting this diagrams in two parts, each with the exter- 
nal vertex p_dxP+ carrying the external frequency lo (see 
Fig. [7]). We should take the product of spectral weights of 
the Green's functions Im g{xi,qi) at the cross section, to 
divide each individual spectral weight by 2 sinh(a;i/2r) 
(in our case of bosons), and to multiply the whole prod- 
uct by 2 sinh(ti;/2T). The sum over all possible cross sec- 
tions, multiplied by the square of the generalized vertex 
part, gives a value of the imaginary part of the diagram. 
See Figures [S] and [HI schematically showing this method. 

Using this method, all contributions can be classified 
by the number and types of Green's functions appearing 
in the cross section. Formally, we have the contribu- 
tions with two propagators in the cross section (2-CS) 
and three propagators in the cross section (3-CS). How- 
ever the cross sections with two propagators are essen- 
tially reduced to the simple loop diagram, with a static 
vertex correction (except for diagrams with ai vertex, 
see below). Given the structure of the external vertex, 
P-dxP+, and kinematic properties for linear spectrum, 
one can verify that the contribution of 2-CS to a does 
not contain cj^^. We may hence exclude this type of 
cross-section from our analysis. 

Each partial 3-CS contribution to the drag resistivity 
can be represented as^ (cf. Eq. ([77)1 ) : 



6ro = nUf^K+K^ do \r 



-jki\ Sjku 



dqidxi 



(90) 



(91) 



S. 



]ki 



g'J {xi, qi)g'l{x2, q2)9i {x3, 93) 



Tsinhl^ 



sinh^ 



sinh ^ 



and the spectral weight g'-{x,q) — ^{q/2)d{qvj ± oj) is 
imaginary part of the retarded Green's function for chi- 
ral density j = R±,L±, Eq. (j2ip . A symmetrization 
prefactor 1 /2 should be added to Sjki when values of two 
indices are coinciding. The total drag resistivity is given 
by summation over all possible j, k, I = R±,L± so that 



E 

jki 



{2itY 

g'j{xi,qi)gk{x2,q2)gi{^ 



dqidq2dxidx2 2 
lljfcil 



Xl - X2-,- 



(92) 
91 - 92) 



2T sinh sinh ^ sinh "^^tg-" 
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R-(L-) 



R-(L-) 



R-(L-) 



FIG. 10: Singular diagrams. 



FIG. 8: Leading contribution to low-T drag coefficient is given 
by diagrams with three propagators in a cross section. 





FIG. 9: First-order in curvature contribution to the vertex 
part in the above diagram. 



We can distinguish between the cross sections with i) 
three different chiral densities, ii) two different chiral den- 
sities and iii) one sort of chiral density. 

I. One sort of chiral density in the cross section 
can appear only on 3-CS of the type and 
g'l_g'l_g'[_- The kinematic restrictions ((9T|) obviously 



lead to the absence of such contribution in ([92|l aX uj ^ Q. 

Having this contribution excluded, we observe that the 
Lo^"^ contribution to a is delivered by setting uj — Q in 
the integrand of (IM)) : upon doing this the calculation is 
simplified. 

II. It can be shown that 3-CS with three different prop- 
agators may contain only two combinations g'^j^g'^+QL- 
and g'fi^gL+9R-- Summing up all the contributions, and 
using the kinematic restrictions in do and Sijk, Eq. (|9ip . 
we obtain that the corresponding static vertex is 



R+L+L- 



27r 



2v+{v 



(93) 

here q refers to L_ boson ; one has Trj^^j^r^ — 
—Tr+l+l-. The contribution to the drag coefficient is 
proportional to 



dqq^\TR,+L+L-\'^{v\ - 



Tsinh 



4T 



sinh 



AT 



'— sinh ^ 



(94) 



this expression corresponds to (jTS]) . However, substitut- 
ing ([88|) into ((93|) . one obtains Tr+l+l- = and the 
absence of this partial contribution to 7'o (cf. next sub- 
section). 

III. The remaining combinations, i.e. 3-CS with two 
different propagators, are : gR-dR-dl-^ gR+9R+9R-^ 
gR+gR+g'L- ^^'^ the combinations which are obtained 
from the above ones by the change R ^ L. Due to obvi- 
ous symmetry, the latter change together with qi ~qi 
leaves the expressions for Stq intact, so it suffices to con- 
sider only the above three CSs and multiply the result 
by two. 



The first cross section, g'^_g'^_g'l_ 
non- vanishing vertex 



corresponds to the 



R-R-L~ 



2ck5W 



27r vl 



(95) 



where the kinematic restrictions were used again. 

Finally, we come to the most problematic term of the 
cross sections g'^_^.g'^_^_g'^_ and g'^^g'^^g'[_- The corre- 
sponding vertices F's in Fig. [51 explicitly read as 



^R+R+R- = aiq3gR+{xi + X2,qi + q2) (96) 
+a2q3gL+ixi + X2, qi + q2) 

-0LiqigR-{x2 + 2:3,92 + 93) 

~'a^qigL-[x2 + xz,q2 + 93), 
Yr.+R+L- = a.iqzgR+{xi + X2,qi + q2) (97) 
+a2qigL+{xi + .T2, qi + 92) 
~oi!,qigR-{x2 + .X3, 92 + 93) 
-0LiqigL-{x2 + xz,q2 + qz)- 

It is seen here, that the term ~ ai is divergent on the 
assumption of ideal bosons, which happens because of 
simultaneous conservation of energy and momentum for 
the linearized spectrum. This divergence in gR+{xi + 
2^2 , 9i +92) would not happen if we would assume a damp- 
ing of boson excitations ; moreover, the kinematic restric- 
tions in the cross section (PT|) lead to prefactor (73 = at 
cj = 0, so that the contribution from ai term is zero 
for any model damping of bosons. In our problem how- 
ever, this "damping" can only be obtained as the result 
of resummation of higher order corrections in . The 
divergent diagrams are depicted in Fig. 1101 

The first diagram here, ^ a\, is the lowest order con- 
tribution to the self-energy part, Eq. pS]). Evidently 
it corresponds to the triple pole {x — v+q)~^ which re- 
sulted in the above divergence. At the same time, as 
discussed above, this is the first singular diagram in a 
series, which should ultimately produce the broadened 
propagator, Eq. (j24p . Importantly, this type of singular 
diagrams appears without the properly fermionic interac- 
tion, as seen in ((88| . Therefore, in the application of the 
formula (|92p . one should omit the most singular contri- 
bution, ~ a\, in |Ffl+fl+fl_p and use the the broadened 
propagator, Eq. ([M]) instead. 

The subleading singular term, ^ ai, is depicted by 
the right diagram in Fig. [TOl We argue here that 
this diagram should vanish in our consideration, and 
provide the following heuristic argument. The linear 
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dispersion law for bosons led to the singularity in ai 
term (1^5)) . but this singularity is removed on the as- 
sumption of finite width of the bosonic pole. In this 
case the Green's function gn^^xi + X2,qi + ^2) in (jM)) . 
(|97p becomes purely imaginary quantity. Considering the 
purely real \Tb,+r+r-\^ = y*R+R+B.^T r+r+b- one sees 
that the linear-in-ai terms disappear. This argument is 
most obvious for the simple Lorentzian form of gR^{x, q) 
and should be generalized for more complicated cases, 
particularly, for gR+{x,q) given by The analytical 

structure of more complicated diag rams at finite T be- 
comes rather involved, see Ref. |23[ and a general proof 
of the cancellation of these vertex corrections should be 
considered elsewhere. Here we make a conjecture that 
the singular terms represented by the second diagram 
in Fig. [TU] vanish exactly for the "dressed" propagators, 
given by p4)) . 



ro - Ul^-fT^, 7 = 71 +72 + 73, (98) 
^ 1 [a3{v+ + v^) ~ a5{v+ - v~)f 

1 alv\ 

with aj given by ([55]) . It is seen that at v+ — V- 
the leading term in 7 is given by term and is oc 
(u-i- — v^)^^, which is different from the above regime 
^ (f+ — u_)"2 non-equal wires, Eq. ([79]) . To resolve 
the seeming discrepancy, we perform the analysis of the 
nearly-equal wires below. 



This conjecture basically amounts to the following con- 
vention : 

i) one should discard the curvature term ai from the 
Hamiltonian, 

ii) one should use the dressed propagators (|24p for 
the modes R+^Lj^, with the "effective mass" 1/m* — 
Q;i/(27r2), and 

iii) the other curvature terms a2 , • ■ ■ , 015 , not leading to 
singularities, can be treated in the regular PT. 
Proceeding this way, we conclude that at higher tem- 
peratures, r > Ti we can use Eq. (175)1 with the change 
Im Xilm X2 —> Im X+Im The finite width of the 
dressed propagators gR+,gL+ leads then to the regime 
ro ~ r^, similar to Eq. 



The interaction-induced corrections^ to the shape of 
propagator (|24p are the next-order effect in curvature, 
and do not influence our calculation of tq below. Indeed, 
at lowest temperatures, T ^ Ti, the peaks in spectral 
weights Im Xi in C3) do not overlap, and it suffices to 
consider the "tail" of one propagator in ([75)1 and the core 
of another one, approximated by (5— function. At higher 
temperatures, T >Ti, propagators in ([75]) significantly 
overlap; this is non-perturbative regime for bosoniza- 
tion, which cannot provide a prefactor in the estimate 
To ~ ■ However, one may argue that all interaction- 
induced details^ in the shape of the core of propagators 
are smeared by higher temperatures and one can safely 
use the free-fermion estimate here. 



At lower temperatures, T < Ti, we obtain the finite 
contribution ro ^ , due to vertices 0:2,..., 0:5, ab- 
sent in the non-interacting situation. A rather long and 
straightforward calculation gives a prefactor in this 
temperature dependence : 



C. Nearly-identical wires 

We consider now the case when the plasmon velocities 
wi_2 in individual wires are nearly coinciding and their 
difference is of order of U12 ■ The main idea of this subsec- 
tion is that the coefficients in the canonical transforma- 
tion diagonalizing the quadratic Luttinger Hamiltonian 
are sensitive to the ratio Ui2/ {vi — v^)- These coefficients, 
in turn, define the values a, for the vertices ([57|) . We can 
calculate these values and use the previous expressions 
((94)) . ((98)) to obtain the drag in this more general case. 

Using dimensionless difference of plasmon velocities 
()67p . it is convenient to introduce the parameter A as 
follows : 



u = A sin A, vi — V2 — w A cos A, 



(99) 



where the first equality is the properly definition of A. 

The general Bogoliubov transformation, diagonalizing 
the quadratic part of the Hamiltonian ([7]), explicitly 
reads as 



with 

^1 - 
^3 = 

S2 = 



diag 
diag 



Si.S2-S3 



1/2 -1/2 1/2 1/2 




(100) 



J2 



Jl ' ^^72 



1/2 1/2 -1/2 -1/2 
,V_ , f + 



^cos(A/2), -sin(A/2), 0, 
sin(A/2), cos(A/2), 0, 
0, 0, cos(A/2), 






'Sin(A/2) 



0, 



0, 



sin(A/2), cos(A/2) 



where v+ and v- are given by (j66p and correspond to the 
new modes 0+ and , respectively. Using this transfor- 
mation, we can determine the values of the coefficients 
ai in (|57)) . 
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Note that the set of the vertices ((87|) was obtained for 
the case A = 7r/2. For other values of A, we have also 
another set of the vertices which is obtained from ([87|) 
by a change i?+ ^ i?_ and L_, explicitly 



TTl 

cur 



3 V 



Rt 



L^Ll) + f{R_Ll 



+a'^{R- + L^)R+L^ 



L-RX) 

(101) 



For our purposes here, it is sufficient to let vji — vj2 
and si ~ S2 in the calculation of aj,a'j . Thus we ne- 
glect the difference in plasmon velocities in (llOOp . At 
the same time, we allow for the intrawire interaction, 
K = vji/si ^ 1. 

We obtain the expressions for vertices in (|57|) as 



0L\ 

a5 



7r2 


3 + 




TO 


2v 








- 1 


TO 


V 




7r2 




+ 3 


TO 










- 1 


TO 


2v 





cos' 



A 
2 
A 
2 



sm 



cos ■ 



2 
A 



cos ■ 



sm — ) sin A, 



sm — I sin A, 
2 ' 



(102) 



a4, 



and the set a'j is obtained from (|102[) by changing aj 
a'j and A — > —A. 

One can verify that despite more complicated struc- 
ture (jlOOp . the drag current vertex is still given by the 
last expression in ((89)) . Further, the existence of the sec- 
ond set of vertices, a'j, leads to the same set of diagrams 
as before, with appropriate change R+ <-> i?_ etc. Im- 
portantly, there are no diagrams of the type Fig. \S\ con- 
taining both a'j and aj. It means that one can use the 
previous expressions (|98p . corroborated by the previously 
vanishing contribution (|94p and to add the contribution 
from the a'^ set, obtained by A — > —A. 

Combining all terms, we find 



(l-i^2)2cos2A + 



Ul^T"^ /cos^A 



4A3 



A2 



0(1) 



sin^ A 



0(1) 



(103) 



This is the central result of this section. Let us discuss 
it in more detail. 

First of all, we see that the term ^ A^^ in (|103p is 
proportional to intrawire interaction strength, (1 — K)'^, 
whereas the term ~ A~'^ is present even for initially non- 
interacting wires. This is not surprizing, because in the 
latter case the factor sin^ A ~ 1/12/ (^i — 1^2)^ means that 
U12 serves as the Luttinger-type interaction {K — 1), this 
term should be obtained for non-equal wires in the order 
of conventional fermionic formalism. 



For identical wires, we have an interesting crossover 
from ?7i2^^ regime at high temperatures to /\Ui2 \ be- 
havior at lower temperatures. Recalling our definitions 
of the drag resistivity R12 = t^tq/ (kpikp2) and of the 
small scale Ti = Ep/S. ~ Ui2kp, we obtain 



Ri 



2^5 



vpElT, ' 



T<Ti^ Ui2kp, 



(104) 



vpsy 



T>Ti, 



which shows a smooth crossover at T ~ Ti. The para- 
metrically large prefactor, l/?7i2, is similar to above Eq. 
(|7ip . This prefactor i) shows non-perturbative charac- 
ter of the obtained result and ii) does not mean actual 
enhancement of the drag effect, due to its overall small 
amplitude ~ T^. 

For the wires with different plasmon velocities we have 
Ti ~ {vi — V2)kp, and the above Eq. (|103p shows the 
following behavior 



Ri 



2^5 



vpElTt 

vlEp' 



T < Ti, (105) 



T > Ti. 



and two regimes do not match at T = Ti. As we saw 
above, regime changes to exponential decrease at T < 
Ti, which means that there is a logarithmically narrow 
crossover region of temperatures, where R12 interpolates 
between and p^ . 



VII. SUMMARY AND CONCLUSIONS 

In this paper we discuss the effects of curvature for 
the Luttinger liquid. The curvature is irrelevant per- 
turbation, in the sense that it does not change the low- 
energy of the system and usually can be simply discarded. 
Bosonization technique regards the curvature as inter- 
action and, at the first glance, a regular perturbation 
theory is possible. However, a certain care must be ex- 
ercised with such PT, in spite of evident correspondence 
between fermionic and bosonization results. We show, 
that the PT in curvature in bosonization produces both 
singular and regular contributions to different quanti- 
ties. The singular contributions arise already for the free 
fermions and correspond to the non-Lorentzian charac- 
ter of fermionic density propagator. We show that these 
singular contributions stem from one vertex of boson in- 
teraction, and we propose to exclude this vertex from PT 
analysis, while simultaneously using the free-fermion ex- 
pression as a "dressed" form of bosonic propagators. For 
the Coulomb drag problem, the dressed form of propa- 
gators becomes important at high temperatures. 

All other bosonic vertices appear due to fermionic in- 
teraction and yield regular contributions in bosonization 
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PT. They are responsible for the low-T regime in the 
Coulomb drag problem. 

Comparing bosonization to fermionic approach at the 
example of Coulomb drag problem, we observe a comple- 
mentary character of these methods. High-temperature 
regime is obtained most easily with fermions, but 
this is essentially non-perturbative regime for bosoniza- 
tion. Conversely, the low-temperature drag regime for 
different wires is obtained relatively easily with bosons, 
whereas the calculation by fermionic approach is cum- 
bersome. We demonstrated that the bosonization is par- 
ticularly useful in case of Coulomb drag between nearly- 
equal wires, when it simultaneously resums several RPA 
channels for fermionic diagrams. 

We did not include spin into our consideration. It also 
requires a separate analysis, because the cubic terms in 
the fermionic densities violate the spin-charge separation 
realized at the level of quadratic action. It is evident that 
spinful electrons in one wire are equivalent to spinless 
fermions in two wires. In our formalism, the Hamiltonian 
is given by with identification (p^ {(p_) for charge 
(spin) density. The above case of nearly identical wires 
corresponds to Zeeman-splitting of the Fermi velocities 
in the magnetic field by value ~ B/kp- 

Finally, we note that in one spatial dimension, the al- 
ternative source for the "ideal" drag was proposed by 
Nazarov and Averini^ They considered the "backscat- 
tering" between the wires, i.e. 2kp Fourier component of 
the interwire interaction t/i2(2fci?) = f^bs, not considered 
here. The effect of renormalization leads to the increase 
IJbs = UbsiEp /Ty~^^ with lowering T, and the drag re- 
sistivity i?i2 ~ U^s diverges in the limit T ^ 0. However, W{x) ~ 
as was noted inii^i, for large interwire distances D, the 
backscattering amplitude is small [/fl ~ ^-^kpO ^^^^ ^^le 
effect of renormalization can compensate this smallness 
only at exponentially small temperatures. The backscat- 
tering drag mechanism is also sensitive to the difference 
of fermionic densities in two wires^ii!^ Particularly, a 
straightforward calculation by Fuchs et ali^ showed the 
exponential suppression of R12 due to backscattering at 
T < VF^kpi — kF2) = In view of an obvious simi- 
larity between and our Ti, it is tempting to regard 
the regime R12 ~ as the only one, surviving at low 



temperatures. The inclusion of the backscattering drag 
mechanism into our discussion is clearly beyond the scope 
of this study. 
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APPENDIX: OPTICAL CONDUCTIVITY 

We calculate optical trans-conductivity according to 
Kubo formula, and have Re aji = iu;~^Im Qji{uj) 
where Qji{ijj) is a retarded current response function. 
Combining different contributions to Q{ijj), depicted in 
Fig. [3] and making analytical continuation to real fre- 
quencies at finite T, we obtain 



, , , N(xi)-N(x2) , , , . 
dqdxidx2Fji ^ ^' \^W{xi)W{x2), 



' uj — xi + X2 + iO 



F12 = F21 = -q'^u\lvl{xi+X2f, 

„2(2 „.2 2\i2 „,2„2\ 



(A.l) 



771 A ( A J. A\j A A A\( , A A\ 

F\\ = <? \^\-'V2q )(X2-V2q )[xiX2+Viq ) 

+ lq^u^vlvl{xl +xl- 2vlq^), (A.2) 



1 



■ {5{x + ei) + 5{x ~ £1) 



2x{e\ - el) 
-5{x + £2) - 5{x - £2)) , 



(A.3) 



and _F22 is obtained from Fn by interchanging vi and V2 
; here N{x) is Planck function. 

We calculate Im Q12 at w 7^ 0. For non-zero frequen- 
cies the coinciding arguments xi ^ X2 — ±£1,2 do not 
contribute and the rest of the integration is simple. Af- 
ter collecting all the terms coming from (jA.3|l we obtain 
the real part of the conductivity Re ct{lj) in the form 
([64)) . with ad given by ((68)) . 
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